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VOLUME INTEGRAL EQUATIONS FOR ELECTROMAGNETIC 
SCATTERING IN TWO DIMENSIONS 

MARTIN COSTABEL, ERIC DARRIGRAND AND HAMDI SAKLY 


Abstract. We study the strongly singular volume integral equation that describes 
the scattering of time-harmonic electromagnetic waves by a penetrable obstacle. 
We consider the case of a cylindrical obstacle and fields invariant along the axis of 
the cylinder, which allows the reduction to two-dimensional problems. With this 
simplification, we can refine the analysis of the essential spectrum of the volume 
integral operator started in a previous paper (M. COSTABEL, E. Darrigrand, 
AND H. Sakly, The essential spectrum of the volume integral operator in electro¬ 
magnetic scattering by a homogeneous body , Comptes Rendus Mathematique, 350 
(2012), pp. 193-197) and obtain results for non-smooth domains that were previ¬ 
ously available only for smooth domains. It turns out that in the TE case, the 
magnetic contrast has no influence on the Fredholm properties of the problem. As 
a byproduct of the choice that exists between a vectorial and a scalar volume inte¬ 
gral equation, we discover new results about the symmetry of the spectrum of the 
double layer boundary integral operator on Lipschitz domains. 


1. Introduction 

The scattering of electromagnetic waves by a penetrable object for a fixed frequency 
is described by the time-harmonic Maxwell equations, valid in the whole space, even 
if the scatterer is confined to a bounded region. For purposes of numerical modeling, 
one prefers a problem posed on a bounded domain, and various methods are known 
to achieve this. One way is the method of volume integral equations (VIEs) that will 
be described in Section [31 This method has been used for a long time in computa¬ 
tional physics O El], but its mathematical analysis is still incomplete. The case 
of coefficients, permittivity e and permeability /i, that are smooth functions on the 
whole space has been studied more completely, because the VIE can then be reduced 
to standard Fredholm integral equations of the second kind, see [SI in HE! and the 
chapter on Lippmann-Schwinger equations in the classical monograph [5j. 

We are rather interested in the case where the coefficients are discontinuous on the 
boundary of the scatterer, as is the case at the interface between different materials. 
The well-posedness of the strongly singular VIEs is then a non-trivial problem. There 
exist negative or complex-valued coefficients e and fi that can nowadays be realized 
with meta-materials, for which the scattering problem and its associated VIE are 
not well posed in the Fredholm sense. For the scattering problem, there has been 
recent progress in finding necessary and sufficient conditions for Fredholmness using 
variational techniques, see [2] for the two-dimensional situation. For the VIEs, the 
question amounts to determining the essential spectrum of the volume integral op¬ 
erators. The knowledge of the essential spectrum is also of basic importance for the 
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analysis of numerical algorithms, for example for the construction of efficient iterative 
solution procedures, see [20J. 

In [2], we analyzed the VIE using a technique of extension to a coupled boundary- 
domain integral equation system, which eventually allows the reduction of the prob¬ 
lem to the better known analysis of boundary integral operators. In the three- 
dimensional case studied in [2], the analysis of the “magnetic” volume integral opera¬ 
tor poses particular difficulties, which is the reason why we could describe its essential 
spectrum only under smoothness assumptions on the boundary and why the analysis 
given in [T51 [IB] is incomplete for this case. 

In this paper we study the two-dimensional case, where it turns out that this 
“magnetic” operator B £ has a very simple structure, see subsection 14.1.21 below, 
in particular Proposition 14.51 This structure implies that the operator, while still 
strongly singular and therefore contributing to the principal part of the VIE, does 
not contribute in a non-trivial way to the essential spectrum. It also suggests a new 
idea for the analysis of this operator in the three-dimensional case that works also 
when the boundary is only Lipschitz, see m Section 2.4] and our forthcoming paper 
on the subject. Thus, while the 2D case has its own interest, for example for the 
study of waveguides or photonic fibres, our main motivation for studying it here is 
that it allows a refined spectral analysis of the volume integral operators. 

Another interesting feature of the two-dimensional electromagnetic problem is that 
one can choose to describe it either as a vector-valued or as a scalar problem. This 
leads to two different VIEs. In both cases, we prove necessary and sufficient condi¬ 
tions on the coefficients e and /i for the volume integral operators to be Fredholm. 
These conditions involve the essential spectrum of certain well-studied scalar bound¬ 
ary integral operators, and at first it looks like the two formulations lead to different 
conditions. They are equivalent, however, because of a symmetry property of the 
essential spectrum of the double layer boundary integral operator that is specific to 
dimension two. We prove this symmetry in Section [5] by relating the boundary in¬ 
tegral equation to a scalar transmission problem involving the coefficient e and for 
which one finds a symmetry with respect to replacing e by A 

2. The two dimensional Maxwell equations 

We consider an electromagnetic scattering problem described by the time-harmonic 
Maxwell equations in M 3 

(2.1) curl E — iu[iH — 0 ; curl H + iueE = J. 

The current density J (which may be identically zero) is related to an incident held 
(E°, H°) by the requirement that the latter satisfies the free-space Maxwell equations 

(2.2) curl E° — iuj^oH 0 = 0 ; curl H° + iui£ 0 E° = J 

with constant permittivity £q and permeability /j,q. In the Maxwell equations (I2.1jh 
the fields are supposed to be locally square integrable, and the equations are satisfied 
in the distributional sense. This means that on a surface of discontinuity of the 
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coefficients /i and e, the fields satisfy transmission conditions 

(2.3) [n x E] = 0; [n x H] — 0; [n • /ii7] = 0; [n • eE] = 0 . 

The scattering problem is completed by the requirement that E — E° and H — H° 
satisfy a radiation condition at infinity. Note that the incident field may contain a 
part generated by the current density J and satisfying the radiation condition, as well 
as a part satisfying the homogeneous Maxwell system but not the radiation condition. 

We assume now that the penetrable obstacle is of cylindrical form 12 x R, where 12 
is a bounded Lipschitz domain in R 2 . In addition, we assume that the fields E and 
H as well as the incident fields do not depend on the variable x 3 . We introduce the 
relative permittivity and permeability such that e = £ r £o, k = HrH o- For simplicity 
we further assume that J has compact support disjoint from the scatterer 12. 

Introducing the wave number k such that k 2 = uj 2 £oHo and suitably renormalizing 
the fields, the resulting system of partial differential equations in R 2 can be written 
as 

d 2 E 3 - i kHi = mi ; d 2 H 3 + i kE x = jj ; 

(2.4) —diE 3 - i kH 2 = m 2 ; -d x H 3 + i kE 2 = j 2 ; 

diE 2 - d 2 Ei - i kH 3 = m 3 ; d\H 2 - d 2 H 3 + i kE 3 = j 3 . 

where 

(2.5) (mi,m 2 ,ro 3 ) T = ik(fi r — 1 )xqH, (ji, j 2 , kY = ik{\ - £ r )xnE + J 

and xn denotes the characteristic function of 12. 

It is well known (and easy to see from (12.4ft ) that this system can be written as 
two uncoupled systems corresponding to the transverse electric (TE) and transverse 
magnetic (TM) polarizations. 

In the TE case, only the transverse components E 3l E 2 of the electric field and the 
longitudinal component H 3 of the magnetic field appear: 

d\E 2 - d 2 E x - i kH 3 = m 3 ; 

(2.6) d 2 H 3 + i kEi = j\ ; 

-diH 3 + i kE 2 = j 2 . 

Note that we understand (12.6j) in the distributional sense, which implies that the 
correct transmission conditions for H 3 and the 2D tangential component of E are 
satisfied on T = <912. 

The TM case is described by the remaining 3 equations of the system (12.4(1 : 

d 2 E 3 — i kH] = mi ; 

(2.7) -<9i E 3 - i kH 2 = m 2 ; 

diH 2 - d 2 Hi + % kE 3 = j 3 . 

We observe the symmetry between (12.6ft and (12.7(1 . which reflects the symmetry 
between E and E[ in the original Maxwell system. It is therefore sufficient to discuss 
one of the two cases in detail, and we choose the TE system (12.6(1 . 
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3. Integral equation formulations for the (TE) case 

There are two ways of reducing the first-order system (12.6j) to second order by 
eliminating one of the two fields, and these two second-order formulations are no 
longer identical, one of them leading to a second order system that can be called 2D 
Maxwell system, and the other one to a scalar Helmholtz equation. Correspondingly, 
we can construct two different volume integral equations, both of them equivalent to 
the TE system (12. bD . but bringing forward different aspects of the volume integral 
operators. 

In the “electric” formulation we have a 2 x 2 system of second order partial dif¬ 
ferential equations, and correspondingly we find a vector-valued VIE. The matrix of 
strongly singular integral operators is the sum of an “electric” operator A ^ depending 
on the permittivity e r and a “magnetic” operator depending on the permeability 
/i r . In the three-dimensional case studied in [9], the electric volume integral operator 
was easier to analyse, using a method of extension to a coupled boundary-domain in¬ 
tegral equations already described in [B]. A much earlier analysis (where this operator 
is called “magnetic”) can be found in |l2j . The magnetic volume integral operator 
is more difficult to analyze in three dimensions, but it turns out that here in two 
dimensions, the operator is a compact perturbation of a projection operator, see 
Proposition 14.51 This simple structure implies that the magnetic contrast does not 
influence the essential spectrum, see Corollary 14.81 

In the “magnetic” formulation, the volume integral operator is a scalar strongly 
singular integral operator, similar to what one gets when acoustic scattering problems 
are treated with the VIE method [Tj. The operator again splits into an electric part 
and a magnetic part, but this time the magnetic operator is weakly singular and 
therefore disappears altogether from the analysis of the essential spectrum. 

For the derivation of the volume integral equation (VIE), we assume in this section 
that the wave number k is a positive real number. Later on, seeing that the operators 
depend analytically on the parameter k, we may consider arbitrary complex values 
for k. 

The TE problem consists of the system of three first order partial differential 
equations (12.61) with compactly supported right hand side, completed by the Silver- 
Miiller radiation condition, which is equivalent to the Sommerfeld radiation condition 

(3.1) lim / \d r u — iku\ 2 ds = 0 

R ~^°° J\x\=R 

for each of the three components u = E\ — E 2 — E ®, H 3 — H® of the scattered held. 
Note that the incident held itself may satisfy the radiation condition if it is generated 
by a current density J that does not vanish identically, but our formulation can also 
handle the case where J is zero and E° does not satisfy the radiation condition, as 
is the case e.g. for an incident plane wave. 

In our scattering problem, the right hand side contains the imposed current density 
and also the induced polarization, which is given by the unknown helds multiplied 
by the inhomogeneity in the coefficients, see equations (12.51) . 
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3.1. The electric formulation: as unknowns. In the following, we use the 

two standard instantiations of the curl operator in two dimensions: 

curlE = d\E 2 — d 2 E\ , curlu = ( d 2 u , — diu) T 
and write (12.61) as 

(3.2) curl E — ikH 3 = m 3 ; curl H 3 + i/cE = j . 

Eliminating H 3 from this system gives a second order system for the unknown E = 

(3.3) curl curl E — /c 2 E = curl m 3 + ikj . 

The construction of volume integral equations is based on the following well known 
result. 

Lemma 3.1. Let gk{x) = jH^\k\x\) be the fundamental solution of the Helmholtz 
equation in R 2 satisfying the outgoing Sommerfeld radiation condition. 

(i) The distribution u satisfies in R 2 the Helmholtz equation 

-(A + k 2 )u = f 

and the radiation condition Mb where f is a distribution with compact support, if 
and only if 

u = g k * f ■ 

If f is an integrable function, then this convolution can be written as an integral: 

u(x) = j 9 k(x - y) f(y) dy . 

(ii) The vector-valued distribution E satisfies in R 2 the system 

curl curl E — fc 2 E = f 

and the radiation condition, where f has compact support, if and only if 

E = fc -2 (V div +k 2 )gk * f . 

Note that the vector part (ii) of this lemma follows from the scalar part (i) and 
the relation 

(V div +k 2 ) (curl curl —k 2 ) = —k 2 (A + k 2 ) 

if one uses the fact that for a function u satisfying the radiation condition also the 
derivatives of u satisfy the radiation condition and that derivations commute with 
convolutions. 

Applying Lemma 13.11 to the system (13.31) and using the system (12 .2[) satisfied by 
the incoming field, we obtain 

E — E° = curl gk * m 3 -- ( V div +k 2 )gk * j 

ik 

where the right hand side is given by (12.51) . but with J = 0. Eliminating H 3 from the 
right hand side using the equation 

curl E — ikp r H 3 = 0 , 
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we finally obtain 

E — E° = curing, * (——-An curl E) - (V div+k 2 )g k * ((1 - e r )y c E). 

Hr 

Thus our VIE takes the form 

(3.4) E - A \E - B u k E = E° 
where we used the abbreviations 

(3.5) rj = 1 - e r , v=l-± 

(3.6) A'lE = -(Vdi v+k 2 )g k * {rj X nE) 

(3.7) B k E = curl g k * {v X ci curl E). 

Strictly speaking, we would not need to write the characteristic function xn here, 
because the coefficient functions 77 and v describing the electric and magnetic contrast 
vanish outside of hi, but this notation allows us to treat g and u simply as numbers 
in the piecewise constant case. 

The VIE (13.4p is valid on all of M 2 . However, as soon as E is known on fi, the 
values of A k E and B k E are known, and then (I3.4j) defines E on the entire space M 2 . 
Thus we will consider the VIE (13.4j) as an integral equation on for the unknown E 
on f], and then use the formula (13.4(1 to extend E to all of M 2 . 

Both integral operators A v k and B k are integro-differential operators that can be 
considered as integral operators with strongly singular kernels, since the second 
derivatives of g k (x — y ) behave like \x — y\~ 2 as a; - y —)■ 0 . 

We use the standard notation for Sobolev spaces H m and for 

H(curl,0) = {E G L 2 {Vt) 2 | curlE G L 2 (Cl)} 

H(div, ft) = {E G L 2 (ft ) 2 | div E G L 2 (ft)} . 

Lemma 3 . 2 . Let 77, v G L°°(ft). Then the operators A v k and B k are bounded linear 
operators in H(curl, H). 

Proof. The convolution with g k is a pseudodifferential operator of order —2 and there¬ 
fore maps any L 2 function with compact support to a function in the Sobolev space 
Hf oc (M. 2 ). For a bounded Lipschitz domain fi, the restriction of this convolution to 
therefore defines a bounded operator from L 2 {Vt) to H 2 (Pt). 

This implies that Af. is bounded from L 2 (fl) 2 to L 2 (Ll) 2 and curl A r k even maps 
L 2 (Q) 2 boundedly into i/ 1 (f2). Hence A k is bounded from L 2 {Ll) 2 to H(curl, Q). 

For B k , we see that it maps H(curl, f2) boundedly to i/ 1 (f2) 2 , hence to H(curl, f2). 

□ 

The main question to be answered here is under what conditions on the coefficients 
77 , v the volume integral operator I — Af k — B k is Fredholm in the space H(curl, fl). 
Equivalently, when is the number 1 not a member of the essential spectrum (see (14.111 ) 
of the operator A k + B k . More generally, the aim is to determine cr ess (H^ + B k ) in 
dependence of the coefficients 77 and u. This is the subject of Section [4] below. 
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3.2. The magnetic formulation: II.> as unknown. We start once more from the 
system (13.2|) and reduce it to a second order equation, this time by eliminating the 
electric field: 

(3.8) — (A + k 2 )H 3 = curl j — ikm 3 . 

The application of Lemma 13.11 leads to 

H 3 - H$ = g k * curl j - ikg k * m 3 = ik curl g k * (j/xqE) + k 2 g k * ((/i r - 1)xqH 3 ) . 
The remaining E on the right hand side is eliminated using 

curl H 3 + ifce r E = J , 


finally leading to 

H 3 -H ° 3 

with 


curl g k * (ax n curl# 3 ) + k 2 g k * (/ 3\nH 3 ) 


(3.9) a = 1 —(3 = fi r - 1. 

Here we used the simplifying assumption that the current density is supported in the 
complement of 12, so that = 0. Otherwise there would appear another term on 
the right hand side. 

We can write this VIE in the form 

(3.10) h 3 - c;h 3 - d*h 3 = Hi 

with the integral operators 

(3.11) C k u — curl g k * (axn curl u) 

(3.12) D^H 3 = k 2 g k * ((3 X nu). 

The natural function space is now simply i/ 1 (f2), and the mapping properties of these 
two integral operators are easy to observe. 


Lemma 3.3. Let a, /3 G L°°(Ll). Then C k is a bounded operator in Lf 1 (f2). The 
operator maps L 2 (Q) boundedly to H 2 (Q), thus it is compact in H l (Vt). 

The question to be answered in the next section is to determine <j ess (C k + D%). 


4. Analysis of the VIE in the TE case 

In order to determine the Fredholm properties of the VIEs (13.41) and (13.101) . we 
study the essential spectrum of the operators A? k + B k and C k + D k using methods 
developed in (SJ EJ 7] and some properties specific to the two-dimensional case. 

As definition of the essential spectrum <r ess (A) of a bounded operator in a Banach 
space we choose 

(4.1) cr ess (A) = {A G C | All — A is not Fredholm} . 

Among the many possible different definitions of the essential spectrum, this one is 
sometimes called the Wolf essential spectrum. It is invariant under compact pertur¬ 
bations. Many of the operators studied in the following, in particular in the case of 
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piecewise constant coefficients, can be reduced to compact perturbations of selfad- 
joint operators in some Hilbert spaces, which implies that if they are Fredholm, they 
are of index zero. 

In this section, we first assume that the coefficients are piecewise constant. This 
means that a, /3, r), v are complex numbers and the volume integral operators depend 
linearly on these numbers: 

K = vK l B" h =vB\ | C? = aCj; Dl = pD\. 

The essential spectrum being invariant with respect to compact perturbations of 
the operators, we can further simplify the VIEs by replacing the fundamental solution 
of the Helmholtz equation by the one of the Laplace equation. 

Lemma 4.1. Let go(x) = — ^ log \x\. Then for any k G C the convolution by Qk — go 
is a pseudodifferential operator of order —4 on M 2 . 

Proof. This can be found in the book HZ), but it is seen most easily by comparing 
the symbols, i.e. the Fourier transforms: 

11 k 2 

w^~w = m 2 -k 2 m 2 

□ 


In particular, the convolution by this difference gk — go maps L 2 (42) boundedly to 
L4 4 (42) for any bounded Lipschitz domain 42. 

We now let the operators A 0 , B 0 , C 0 , D 0 be defined like Af B\, Cf Df but for 
k = 0. More precisely 

(4.2) H 0 E = -(V div)(7o * (xnE) 

(4.3) B 0 E = curl 5 -o * (xn curl E) 

(4.4) C 0 u = curl g 0 * (x.n curl u) 

(4.5) D 0 u= 0. 


Corollary 4.2. The operators 

A\ — A 0 , B\ — Bo : H(curl, 42) —> H(curl, 42) and 

Cl ~ C 0 , D\ — D q \ iL 1 (42) —> iL 1 (42) 


are compact. 


Our task is now reduced to finding the essential spectrum of the operators gAo+isBo 
and aCo + /3D 0 . 


4.1. Analysis of the electric formulation. We want to determine the essential 
spectrum of the operator gA 0 + iyB 0 in H(curl, 42) for any 77 , v e C. In order to 
reduce this problem to a better-known one, we use the technique described in jSJ[9] of 
extension by integration by parts. This transforms the strongly singular VIE into an 
equivalent system coupling a weakly singular VIE with a boundary integral equation. 

We use the definitions (l4.2jMI4.3jh where the derivatives and the convolution com¬ 
mute as operations on functions or distributions on M 2 . The multiplication by the 
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characteristic function is understood as the operator of extension of a function 
defined on fl by zero outside of fl. Derivatives do not commute with y^, but we have 
the formulas of integration by parts or weak form of Green’s formulas 

(4.6) div y^E = xn div E — y'n • E for E G H(div, O) 

(4.7) curl y^E = yn curl E — y'n x E for E G H(curl, fl). 

Here the embedding y' : —> 77“^ (M 2 ) is the adjoint of the trace operator, 

and 

EiAn-E: H(div, D) H~* (r) and 
E4 nxE: H(curl, fi) -A- 77(r) 
are the normal, respectively tangential, boundary trace mappings. 

4.1.1. The electric operator. For E G H(div,h2) we can now write 

(4.8) A 0 E = - V^FdivE + VSSn-E. 

Here jY is the Newton potential 

YYu{x) = / go(x — y) u(y) dy 
Jci 

and SS is the harmonic single layer potential 

SSv(x)= I g 0 (x-y)v(y)ds(y). 

Using the fact that the Newton potential is a right inverse of —A and that single 
layer potentials are harmonic on the complement of T, we see that (14. 8 p implies for 
any E G H(div, ft) 

(4.9) divHoE = divE. 

We now split L 2 (fl) 2 into an orthogonal sum 

(4.10) L 2 (0) 2 = V77 0 1 (Q) ©H(divO,D) 

where the second summand is defined by this orthogonality. Later on in Section Owe 
will also need the orthogonal decomposition 

(4.11) L 2 (D) 2 = V 77 1 (0) © H 0 (curl 0, Q). 

We observe that in (U.lOp both summands are invariant subspaces of A 0 : On (p G 
/7q( 0), V commutes with yn, so that 

H 0 V p = - V div g 0 * (y n V p) 

= - V div V 5-o * {xnp) = V ip . 


That is, A 0 acts as the identity on the first summand in (14.1011 . According to (14.9p . 
Ao maps the second summand into itself. 
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Therefore we only need to study the spectral properties of Aq on H(div0,fl). 
Taking the normal trace of (j4.8j) and using the jump relations of the single layer 
potential gives 

(4.12) n • A 0 E = n • V SSn E = ±n • E + K' n • E 

with the integral operator of the normal derivative of the single layer potential on T 

(4.13) K'u(x) = I d n{x) g 0 (x - y)u(y) ds{y) {x G T). 

The extension method then consists in considering n • E as an independent unknown 
and writing the equation for the resolvent of A 0 

AE — A 0 E = F 

as the system that we obtained for E G H(div0,fl) and v = n • E G 77“^ (T), where 

w = n • F: 

f AE — V SSv \ _ (F\ 
yAu - (|I + K')vJ ~ \ w) 

Due to its triangular nature, this latter system is easy to analyze: The terms on the 
diagonal are All and AI — (|I + K') and their invertibility or Fredholmness determines 
the corresponding property of A 0 . 

As a tool, we need the essential spectrum of the operator d n SS = + K' in the 

space 77~^(T). We write 

(4.14) Z = a ess {±I+K'). 

The following result is known da hue] 

Lemma 4.3. Let T be the boundary of a bounded Lipschitz domain fl C R 2 . Then 

h + K' : 77“^(T) —>■ H^(T) 

is a selfadjoint contraction with respect to a certain scalar product in H~^(T). Its 
essential spectrum E is a compact subset of the open interval ( 0 , 1 ). 

7/T is smooth (of class C 1+a , a > 0), then K' is compact, hence £ = {|}. 

7/T is piecewise smooth with corner angles ujj C (0,27 r), then £ is an interval con¬ 
taining namely the convex hull of the numbers y- and 1 — y. 

As a result of the extension method, we obtain the essential spectrum of Aq. 

Theorem 4.4. Let Vt C ! 2 be a bounded Lipschitz domain. The essential spectrum 
of the operator Aq in L 2 (D ) 2 is given by 

CTess(Ao) = {0} U £ U {1} 

where £ is defined in (j4.14j) . The operator A 0 acts as a bounded operator in each one 
of the spaces H(div, fl) and H(curl, il). In both spaces, the essential spectrum is the 
same as in L 2 (fl) 2 . 
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Proof. A quick way to prove that the extension method described above does indeed 
lead to a characterization of the essential spectrum of A$ is to use the well-known 
“recombination lemma”, namely that for two linear operators 

S : Y -A X and T : X -A Y 

the two products ST and TS are spectrally equivalent outside zero. This means that 
for A 7 ^ 0, the two operators XIx — ST and Aly — TS are simultaneously Fredholm or 
not and have the same kernel and cokernel dimensions. For a proof, see 0 Lemma 
4.1] or [TS], p. 38]. 

In our case, we factorize A 0 in H(div0, h2) as A 0 = ST with 
S : H~^ (r) —>• H(div 0, Q); Sv = V SSv , 

T : H(divO,U) ^hH(r); TE = n • E . 

Thus A 0 is spectrally equivalent outside zero to 

TS = n-V SS = h + K' on H~*( r). 

That A = 0 is in the essential spectrum follows from (14.8[i directly: Aq = 0 on the 
space 

H 0 (divO, U) = {E G H(div, U) | divE = 0, n • E = 0} . 

Thus on H(div0,h2), the essential spectrum is {0} U E, and together with the de¬ 
composition (I4.10|) and A 0 |v^i(n) = ^ we obtain cr ess (A 0 ) = {0} US U {1}. The 
essential spectrum is the same in H(div, fl) and in H(curl, O), because both spaces 
are invariant under Aq, and the operator S defined above maps H~^(T) in fact into 
H(curl, Q) fi H(div0, Q). 

□ 

4.1.2. The magnetic operator. We can proceed along the same lines for the operator 
Bq defined in (I4.3|h Using the integration by parts formula (14.7[) . we get for E G 
H(curl, f1) 

BoE = curl curl g 0 * (xoE) + curlS'S'n x E 
= E - A 0 E + curl SSn x E 

If E G H(curl, fl) D H(div, 0), we can further reduce this using (14.81) to 

(4.15) B 0 E = E + V div E — V • E + curl SSn x E . 

If we want to use the extension technique like we did for A 0 , we can reduce B 0 to 
an equivalent boundary integral equation, but we now need both the normal and 
tangential traces on the boundary. For this we use the orthogonal sum 

(4.16) H(curl,Q) = VF^U) ® X with X = H(curl, Q) n H(div0,U). 

Note that this decomposition is orthogonal both with respect to the L 2 scalar product 
and the H(curl) scalar product. 
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On V Hq(Q) we have Bo = 0 by definition. On X, we use (I4.15|) and the recombi¬ 
nation lemma and factorize I — Bq — ST, where S and T are now defined as 


S : H— 2 (T) 2 X; S[ = V SSv * curl SSw , 


T : X —>■ (T) 2 ; T E = 


n • E 
n x E 


The recombined operator TS is a 2 x 2 matrix of boundary integral operators involving 
the normal and tangential derivatives of the single layer potential: 

K' = n • V SS = -n x curl SS 


§ 1 - 


and 


(4.17) 


TS = 


d t SS = n x VSS = n • curl SS 

Thus I — Bq on X will be spectrally equivalent outside zero to the matrix operator 
on tH(T ) 2 

'jl + K' -d t SS' 

'dtSS |I + 7 

Even for a smooth boundary T, the essential spectrum of this matrix of operators is 
not obvious, because whereas K' is then compact, the strongly singular integral oper¬ 
ator dtSS is not compact. But for smooth T we can use symbols of pseudo differential 
operators to find the essential spectrum of the operator matrix, and for piecewise 
smooth T we can use Mellin transformation to the same end, but we have not seen 
an obvious way to find the essential spectrum of TS for a general bounded Lipschitz 
boundary T. 

We will instead turn around the equivalence between I — B 0 and TS by finding a 
different way to analyze B 0 and thence deduce, as a corollary, the spectrum of TS in 

(SUD- 

We go back to the definition (14.31) of B 0 and use curl curl = —A to obtain 
curl Bq E = xn curl E. On the space H(curl, fl), we have therefore 

(4.18) curl B 0 E = curl E . 

This implies 

Bq E = curl go * (xn curl E) = curl g 0 * (xa curl B 0 E) = B 2 E . 

Proposition 4.5. The operator Bo on H(curl, fi) is a projection operator. 

Both the kernel of Bo, which is H(curl 0, Q), and the image of Bq, which is the 
orthogonal complement of H(curl 0, Q), hence contains curlif^fl) H H(curl, O), are 
of infinite dimension. Therefore the essential spectrum is equal to the spectrum: 


a ess (B 0 ) = {0,1} . 

As a corollary, we get the same essential spectrum for the matrix operator TS in 
(14.171) . now valid for any bounded Lipschitz T. But by a careful inspection of the 
recombination argument we can do more: We have B 0 S = 0 and therefore 

(TS) 2 = T(ST)S = T(I - B 0 )S = TS. 
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Corollary 4.6. For any bounded Lipschitz domain O C R 2 with boundary T ; the 
matrix TS of boundary integral operators in (I4.17j) is a projection in H~^(jT) 2 . 

This result is equivalent to the following two interesting relations for the tangential 
derivative of the single layer potential (Hilbert transform) on a Lipschitz boundary: 

(4.19) ( d t SS ) 2 = \1- (A '') 2 and d t SS K' + K'd t SS = 0 . 

4.1.3. The electromagnetic operator. We combine now the analysis for the operators 
Ho and Bo and find the essential spectrum of r/Ao + vBq. As before, we use the 
orthogonal splitting (I4.16j) into electrostatic fields and electrodynamic waves. On the 
first summand V Hq(Q), we have r]A 0 + iyB 0 = rjl. Hence 77 G a ess (r]A 0 + uB 0 ), and we 
only need to consider the operator in the second summand X. We use the extension 
method, and after seeing how this worked for the operator Bq, we do not use the 
extension by the normal and tangential traces on the boundary, but by curlE in fl 
and n • E on T. Thus we write for E G X 

(4.20) rjA 0 + uB 0 — r] V SSn ■ E + v curl g 0 * (xn curl E) 

We see that we can factorize r]A 0 + vB 0 = ST, where S and T are now defined as 

s(l) = r]VSSv + ucurlFKF, 

TE = 

The recombined operator TS now acts on Fd~^(Y) x L 2 (Q) and is given by 

ZiA /7711 • V SSv + un ■ curl xVF\ + K')v + vn ■ curExFF 

\ F J y v curl curl ,JFF ) y uF 

The essential spectrum of the triangular matrix 

tq _ AKP + K') on ■ curlt/L 
1 b ~ \ “ 0 vl 

is obvious: a ess (TS ) = r/EU {u}. To this we might have to add the point {0}. 

Let us show that 0 is not in the essential spectrum if rj ^ 0 and v ^ 0. We know 
that 0 is then not in the essential spectrum of TS because 0 ^ £. Thus TS is a 
Fredholm operator. One can see that T is also a Fredholm operator, because the 
boundary value problem 

div E = 0 and curl E = F in Q 
n • E = v on T 

has, for any F G L 2 (fl), v G a solution E G X, as soon as f r v ds = 0 for 

all connected components Tj of T. The solution is unique up to a finite-dimensional 
kernel. For a detailed proof of the corresponding result in M 3 , see [T, Prop. 3.14]. 
Now if TS and T are Fredholm, then S is Fredholm, too. Hence ST is a Fredholm 
operator, and 0 jL cr ess (ST). 

We have proved the main result of this section 


n • E 

curl E 


S : frs(r) x A 2 (H) ^ X; 
T : X — >■ H~^ (r) x A 2 (H); 
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Theorem 4.7. Let C R 2 be a bounded Lipschitz domain. For any rj, v e C, the 
essential spectrum of the operator ijAq + uBq in the space H(curl, ff) is given by 

a ess (r)A 0 + uB 0 ) = {rj} U 77 S U {v} . 

Corollary 4.8. Let Ft C R 2 be a bounded Lipschitz domain. Let the coefficients 
e r and p r be complex constants. Then the integral operator of the VIE (15.41) is a 
Fredholm operator in H(curl, fi) if and only if 

e r 7 ^ 0 , n r 7 ^ 0 and e r 7 ^ 1 - for all a G £ . 

a 

If the boundary T is smooth, then the condition is 

e r ^ { 0 , — 1 } and fi r 7 ^ 0 . 

Proof. The operator I — A\ — Bf is Fredholm if and only if 1 ^ a ess (A v k — Bf) = 
o’ ess (r]A 0 + vBq). According to the theorem, this is equivalent to the conditions 

1 7 - r), 1 ^ g 1 7 -pa for all a 6 £ . 

With the definitions (15.51) . this is easily seen to be equivalent to the stated conditions 
on e r and p r . For smooth F, we use the fact that £ = {|}. □ 

4.2. Analysis of the magnetic formulation. In this paragraph, we study Fred¬ 
holm properties of the VIE (15.101) for the magnetic held H 3 in i/ 1 (f2). Since the 
magnetic formulation describes the same TE transmission problem (12, 6 [) as the elec¬ 
tric formulation studied in the previous paragraph, we expect equivalent results to 
those given in Theorem 14.71 and Corollary 14.81 

Whereas previously we had a two-dimensional Maxwell system that in some re¬ 
spects closely resembles the full three-dimensional Maxwell system, we have now, 
however, a scalar Helmholtz problem that is specihc to the two-dimensional situa¬ 
tion. For the analysis of the scalar VIE we use similar techniques as before, namely 
the partial integration and extension method, but we can take advantage of the sim¬ 
pler structure and allow more general coefficients. In the following, we assume 

(4.21) aeC(ii); [5 e L°°(Fl). 

Note that according to the definition (15.9|) of a, this already includes the condition 
e r 7 ^ 0 on fh Recall that in the absence of further specification, we assume that 
C M 2 is a bounded Lipschitz domain. 

According to Lemma 15.51 the operator D k is compact. This implies that it, and 
with it the coefficient j3, hence the magnetic permeability fi, will disappear from the 
further discussion. It has no influence on the Fredholm properties of the TE problem. 
Note that also in the previous paragraph the only condition on /i we encountered was 
Hr 7 ^ 0. This was simply a necessary prerequisite for the definition of the coefficient 
v in the electric VIE (15.41) . The condition v 7 ^ 1 that is required according to 
Theorem 14.71 is always satisfied if v is defined from ji r as in (15.5|) . 

Theorem 4.9. Let C R 2 be a bounded Lipschitz domain. Let the coefficients a 
and /3 be defined by 



/ ^ hr 1 












VOLUME INTEGRAL EQUATIONS FOR 2D ELECTROMAGNETISM 


15 


and satisfy the conditions (14.211) . Then the operator I — C k — of the VIE (13.10)) 
is Fredholm in if and only if the boundary integral operator 

(1 - f )I - aK : H 1 / 2 (T) —>• H 1 / 2 {T) 

is Fredholm, where K is the boundary integral operator of the double layer potential, 
see (I4.23P below. 

A sufficient condition for this is that for all x e T and a E £ (defined in ()4.14j ) ) 

erO) ~f~ ~~~r ■ 

a — 1 

If e r is constant, then this condition is also necessary. 

If T is smooth, the necessary and sufficient condition is e r (x) ^ —1 for all x G T. 


Proof. Our task is reduced to the determination of the Fredholm properties of the 
operator I — C k defined in (13.111) : 

C k u = curl g k * (axci curl u ). 

We use the partial integration formula (dual to (14.7)7 

(4.22) curl xnu — xn curl u + y'tyu for u E H 1 (fl), 

where t is the unit tangent vector and 7 u is the trace of u on T. With the product rule 
a curlu = curl(cm) — u curia and the identity curl g k * curl = —A g k * = I + k 2 g k *, 
we can write 

C k u = curl g k * (a curl xn u ) — curl g k * (y'tyaw) 

= curl t/ fc * (curl(axnw)) - curl g k * (xciu curia)) - curlg fc * (y'tyau) 

= au + Cpa + @h(l m f) • 

Here C\ is the weakly singular integral operator given by 


Ciit (a;) = k 2 


g k (x - y)a(y)u(y) dy - curl 


g k (x - y)u(y ) curla(y) dy . 


J ct Jci 

It maps L 2 {kl) boundedly to H 1 ( fl) and is therefore a compact operator in 7/ 1 (0). 
The operator 3> k is the Helmholtz double layer potential 


- curl g*, * (y't jau)(x) = ^( 7 au){x) = j d n ^g k (x - y)a{y)u(y) dy (x e H) . 
Thus 

(I — C k )u = (1 — a)u — C\U — @ k ('yau ). 

Modulo compact operators, we are left with the operator 

u 1 —y (1 — a)u — 3>o(x au ) = j-(u — e r ^ 0 ( 7 cm)) , 


where is the harmonic double layer potential. Since by assumption e r 7 ^ 0 on 0, we 
only need to consider the operator I — e r @o(y a ’)- By the recombination lemma, the 
Fredholm properties of this operator on H 1 (fi) are the same as those of the operator 

I — 7 (ae r ^o) on H 2 (r). 
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Here we used the factorization e r S>o{^a-) = ST with 

S : H*(T)-> H\n)-, Sv = e r %v, 

T : H\Sl) —> H^(T); Tu = ^{au). 

We then use the jump relation 7^0 = — + K for the double layer potential with 

(4.23) Kv{x) = I d n ( y) go(x - y)v(y) ds(y) (iGT). 

Thus 

I — 7(ae r ^ 0 ) = I + 1 — ae r K = e r ((4- + f )I — aK) 

Now 4- + | = 1 — and we see that we have to determine the Fredholm properties 
of the boundary integral operator 

G := (1 - f )I - clK = I - a(§I + K) 

as claimed. Since K is the adjoint operator of K', we have 

a ess (kj + K) = a ess {±I + K') = T. 

If a is constant on T, we see that G is Fredholm if and only if aa 7 ^ 1 for all a G £. 
This is equivalent to the stated condition e r 7 ^ - JL j. If a is not constant on F, we can 
use a standard localization procedure based on the fact that the commutator between 
the singular integral operator K and the operator of multiplication with a smooth 
function is compact in H^{T). If the condition e r (x) 7 ^ -S— for all cr 6 S is satished, 
one can construct a local regularizer using an inverse modulo compact operators for 
the operator with e r frozen in the point x. These local regularizers are then patched 
together using a partition of unity, giving a global regularizer that shows that G is 
Fredholm. We omit the technical details here, see m Section 3]. Finally, if T is 
smooth, then K' is compact, therefore G is Fredholm if and only if (1 — |)I is. This 
is the case if and only if a(x) 7 ^ 2 for all ifT, which is equivalent to e r (x) 7 ^ — 1 . □ 

5. Symmetry of the spectrum 

The condition on e r in Theorem 14.91 e r 7 ^ is not the same as the one given in 
Corollary 14.81 e r 7 ^ 1 — 4. They coincide only if 

uGE 1 - cr G E. 

This symmetry of the essential spectrum of the operator + K' with respect to 
the point \ is a non-trivial result that may be new. It can in fact be seen as a 
consequence of the fact that both VIEs considered in Corollary 14.81 via the electric 
held formulation and in Theorem 14.91 via the magnetic held formulation are equivalent 
to the same Maxwell transmission problem. Since this amounts to a very indirect 
way of proving the result, we shall give an independent proof below that does not use 
volume integral equations. Note that the result is known for smooth domains where 
£ = {|} and for piecewise smooth domains where Mellin analysis shows that £ is an 
interval with midpoint |, but we show it for arbitrary bounded Lipschitz domains in 
M 2 . It is probably not true in higher dimensions. 
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Theorem 5.1. Let 12 C M 2 be a bounded domain with a Lipschitz boundary T. Let 
K be the boundary integral operator of the harmonic double layer potential acting on 
(r). Then 

(5.1) cr ess (K) = cr ess (■- K) . 

This symmetry implies that the operators + K and ~I — K in H^(T) and the 
operators + K' and — K' in all have the same essential spectrum E, 

and the same holds for the Helmholtz double layer potential. 

The proof uses the equivalence between the boundary integral equation and a scalar 
transmission problem together with the following known symmetry result for scalar 
elliptic boundary value problems. 

Lemma 5.2. Let 12 C R 2 be a bounded Lipschitz domain and let the complex-valued 
coefficient function e G L°°(12) satisfy 1 G L°°(12). Then 

(5.2) dive V : H] (Q) -a H~\Q) 
is Fredholm if and only if 

(5.3) divl V : H\n) -A (H 1 ^))' 
is Fredholm. 


This lemma is due to |2J Theorem 4.6], where it was proved under the hypotheses 
that 12 is simply connected and that e is real-valued. We give an independent proof 
that does not need these additional assumptions. Note that, as A : F1 q( 12) —>• if _ 1 (l2) 
is an isomorphism, the operator in (1 5.2 1) is Fredholm if and only if the mapping u i—)■ in 
defined by 

(5.4) / e'Vu-’Vvdx= / Vw-Vucfe for all v G H ( ] (Q) 

Jvl Jn 

is Fredholm in Hq (Q). Similarly, the operator in (15.311 is Fredholm if and only if the 
mapping u i—)• w defined by 


(5.5) 



- V u ■ V v dx 


is Fredholm in // 1 (12). 


V w ■ V v dx 


for all v G H 1 (12) 


Proof. We use a simple linear algebra lemma whose elementary proof we leave to the 
reader. 


Lemma 5.3. Let X, Y, Z be vector spaces and S : Y —>■ Z, T : X —>• Y be linear 
operators such that S is surjective and T is injective. Then T induces an isomor¬ 
phism between ker(ST) and ker S flimT, and S induces an isomorphism between the 
quotient spaces T^kerS 1 + imT) and Z/ im(ST). In particular, ST : X —>■ Z is an 
isomorphism if and only ifY — ker S © im T, and ST is a Fredholm operator if and 
only if 

dim(ker S D im T) < oo and codim(ker S + im T) < oo . 
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We apply this lemma to the situation 

div : L 2 (flf ->■ H~\n) , e V : H%( 0) —>■ L 2 (D) 2 
and obtain that dive V : Hq(Q) —> if _1 (n) is Fredholm if and only if 

dim(H(divO, fl) D eV < oo 

and codim(H(div 0, fl) + eV Hq(Q)) < oo in L 2 (fX) 2 . 

Using the well known fact that there is a finite dimensional space M’q such that 
H(div0,h2) = curliL^Q) + we see that this condition is equivalent to the fact 
that 
(5.6) 

dim(curlif 1 (0) D eV H$(&,)) < oo and codim(curli7 1 (0) + eV Hq(Q)) < oo . 

Dividing by e, we see that the latter condition is equivalent to 

dim(i curl fl V Hq(Q,)) < oo and codim curl + VifpfO)) < oo . 

Using the definition fl4.HH . it is also well known that there exists a finite dimensional 
space M -_( such that H 0 (curlO,n) = V and therefore we obtain the 

equivalent condition 

dim(i curl if 1 (Q)fiH 0 (curlO, O)) < oo and codim(i curlif 1 (O)+H 0 (curl0, fl)) < oo . 
According to Lemma 15.31 applied to the situation 

curl : L 2 (Cl) 2 -A (if 1 (Q)) / , \ curl : H l {Fl) —>■ L 2 (D) 2 

the latter condition is equivalent to the Fredholmness of the operator 

curl i curl : H\Q) -A . 

We finally note that in two dimensions, 

curl - curl = — div - V . 

e e 

Therefore the last condition expresses the fact that the operator in (15.3j) is Fredholm. 

□ 

In order to apply this result to our boundary integral operator, we consider the 
following situation: UTU fi + , where the interior interface F is Lipschitz and 

satisfies T = dfi~ = \ dil. The coefficient function e is piecewise constant, e = 1 

in fD", e = e r e C in fl _ . Define K' as the boundary integral operator of the normal 
derivative of the single layer potential, adjoint of the double layer potential, acting 
in the space H~z(T), see (j4.13j) . In this situation, we have the following result. 

Lemma 5.4. The operator diveV, acting either from H ( \ (Q) to iL _1 (0) or from 
H\Q) to , is Fredholm if and only if the integral operator 

-^2—^ + (u _ 1 )-^ / 

is a Fredholm operator in H~^(T). 
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Proof. We consider the first case (Dirichlet condition) of dive V defined on Hq(Q). 
The second case (Neumann condition) is treated very similarly, the space Hq(Q) being 
replaced by Hl(Q), the space of if 1 (fl) functions of mean value zero. We use the 
variational formulation (15.4p for the operator dive V on Hq(Q). We split the space 
according to 

H^(n) = Hi r (fi) © J^ r (0) , 

where 

Hlr(n) = {ue Hl(n) | w,r = 0} = ^(fT) 0 H^(n + ) 

and the second term is defined as the orthogonal complement of Hl r (Q) in 

Hq(Q). Orthogonality is understood in the sense of the Hi scalar product appearing 
in (I5.4[) . 

The subspace Hl r (Q) is invariant under multiplication by e, and it is also invariant 
under the operator A -1 diveV defined by the variational formulation (15.41) . This 
operator is invertible on this subspace as soon as e r ^ 0. Indeed, if we decompose 
orthogonally 

u = U- + u + + Uq , w = W- + w + + Wo 
with u±,w± G HK^) and Uo,Wo G M'l r (0), then we can partially solve (15.4|) : 

u + = w + , = -u?_ . 

This shows that our operator div e V is Fredholm if and only if in the restriction of 
the variational formulation (15.4|) to 

(5.7) / e'Vuo-'Vvdx= / V wq ■ V v dx for all v G J£^ r (Q) 

Jn Jn 

the mapping uq K > Wo is a Fredholm operator in From the dehnition of 

as orthogonal complement we hnd the characterization of this space 

J#J r (Q) = {u 0 G Hl r (Q) | Au 0 = 0 in fi" U D +}. 

The standard representation formula shows that the functions in this space are rep¬ 
resented as single layer potentials 

u 0 (x) = SS D (j)(x) = j G D {x,y)<f>(y)ds{y). 

Here Go is the Green function of the domain fl, that is the fundamental solution 
of —A satisfying homogeneous Dirichlet boundary conditions on the outer boundary 
<9f 1. It is a classical result that Gn(x,y) = go(x — y) + h D (x,y ) with a C ,00 (D x D) 
function hp. The jump relations for the single layer potential SSp are then the 
same as for the single layer potential SS , in particular for the one-sided traces of the 
normal derivatives df and d~ on V 

d%SS D <i> = Ti<i> + K' D <l> 

where K' D — K' is a compact operator on H~^(T), due to the smoothness of h D . 
The representation of Uq by a single layer potential defines an isomorphism between 
H~z{T) and with inverse given by the jump 0 = d~u 0 — dfuo- 
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We can now use Green’s formula in f2 + and for (j5.7j) and see that this variational 
formulation is equivalent to the transmission condition 

e r d~u 0 -d+u 0 = d~w 0 -d+w 0 . 

Representing uq = SSd4> and wq = SSd^> and using the jump relations, we find the 
equivalent boundary integral equation 

(e r + 1)|0 + (e r - 1 )K' D (j) = 0 . 

We have thus shown that the operator (e r + l)|l + (e r — 1 )K' D is Fredholm in H5 (r) 
if and only if dive V is Fredholm from to its dual space. 

This proves the “Dirichlet” version of the lemma. In order to show the “Neumann” 
version, one can repeat the same steps, replacing 7/^(0) by Hl(Q) and the Green 
function Gd by the Neumann function Gn, that is the fundamental solution satisfying 
homogeneous Neumann boundary conditions on dVt. One obtains a boundary integral 
equation involving a “Neumann” version K' N of K' , which will also differ by a compact 
operator from K' and thus have the same essential spectrum as either K' D or K'. □ 

Combining Lemmas 15.21 and 15.41 it is now easy to finish the proof of Theorem 15.11 
We find that for any e r G C \ {0}, the operator 

(e r + l)yl + (e r — 1)A 

is Fredholm in H~^(T) if and only if the operator 

is Fredholm. Setting e r = we have proved that for any A ^ ±|, AI + K' 

is Fredholm if and only if AI — K' is Fredholm, which concludes the proof of the 
theorem. 
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